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Abstract 

We consider a two type (red and blue or R and B) particle popula- 
tion that evolves on the d-dimensional lattice according to some reaction- 
diffusion process R + B — > 2R and starts with a single red particle and 
a density p of blue particles. For two classes of models we give an upper 
bound on the propagation velocity of the red particles front with explicit 
dependence on p. 

In the first class of models red and blue particles respectively evolve 
with a diffusion constant Dr = 1 and a possibly time dependent jump 
rate Db > - more generally blue particles follow some independent 
bistochastic process and this also includes long range random walks with 
drift and various deterministic processes. We then get in all dimensions 
an upper bound of order max(p, ^fp) that depends only on p and d and 
not on the specific process followed by blue particles, in particular that 
does not depend on Db- We argue that for d > 2 or p > 1 this bound 
can be optimal (in p), while for the simplest case with d — 1 and p < 1 
known as the frog model, we give a better bound of order p. 

In the second class of models particles evolve with exclusion and possi- 
bly attraction inside a large two-dimensional box with periodic boundary 
conditions according to Kawasaki dynamics (that turns into simple exclu- 
sion when the attraction is set to zero.) In a low density regime we then 
get an upper bound of order y/p. This proves a long-range decorrelation 
of dynamical events in this low density regime. 



1 Models and results 



1.1 A diffusion-reaction model 

In [6] Kesten an Sidoravicius considered the following Markov process. A count- 
able number of red and blue particles perform independent continuous time sim- 
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pie random walks on the d-dimcnsional lattice Z d . Red particles jump at rate 
Dr and blue particles jump at rate Db- When a blue particle jumps on a site oc- 
cupied by a red particle, the blue particle turns red. When a red particle jumps 
on a site occupied by blue particles these turn red. Thinking respectively at the 
red and blue particles as individuals who have heard about a certain rumor and 
are ignorant of it - or as individuals who have or have not a certain contagious 
disease - this Markov process provides a model of rumor propagation - or epi- 
demic diffusion - inside a moving population. This is also a reaction-diffusion 
dynamics of the kind R+ B — > 2R that can model a combustion process. 

We define a each time t > a red zone lZ(t) that is the set of sites Z d that 
have been reached by some red particle at some time s £ [0;(]. At any time 
t > all the red particles stand in the red zone, but some blue particles can 
stand in the red zone and the red zone can contain empty sites. The red zone is 
the set of sites reached by the rumor or the set of burnt sites according to one 
or another interpretation of the process. 

Let us assume that the initial configuration was built in the following way. 
We put independently in each site z £ 7L d a random number of blue particles 
according to Poisson variables of mean p > 0, then at time t = we choose 
one particle according to some probabilistic or deterministic rule, we turn it red 
and we turn red the possible other particles that stood in the same site. Then, 
denoting by B(z 1 r) the Euclidean ball of center z and radius r and making a 
change of origin to have 71(0) — {0}, Kesten and Sidoravicius proved [5J: 

Theorem [Kesten-Sidoravicius]: If Db = Dr > there are two positive 
and finite constants C\ < C 2 such that with probability 1 



will hold for all t larger than some finite random time Tq . 

If Dr > there is a finite constant C 2 such that with probability 1 



will hold for all t larger than some finite random time To . 

Remarks: i) Actually they did not introduced any change of origin. The 
analogous result without change of origin is an equivalent statement, but our 
change of origin will serve us later. 

ii) They proved the theorem in a slightly more general situation: when the 
initial configuration is obtained by adding any finite number of red particles in 
a finite set set of sites to a Poissonian distribution of blue particles. However it 
is easy to see that the same result in this more general case is equivalent to the 
previous theorem. For the sake of simplicity we will restrict ourselves to discuss 
processes built like above, starting with a single blue particle that turns red. 

iii) The inclusion p.2[) gives a "ballistic upper bound" on lZ(t). The "ballistic 
lower bound" expressed in (jl.ip is much harder to prove and was obtain only in 
the special cases Db = Dr > [6J and Ab=0([2], [3], 0])- But it is believed 
that such a bound holds in the general case Dr > (see [5J). 



B(Q,dt) C TZ(t) c B(0,C 2 t) 



(1.1) 



TZ(t) c B(0,C 2 t) 



(1.2) 
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iv) On the basis of (jl.ip . that is of a ballistic upper and lower bound on TZ(t), 
Kesten and Sidoravicius proved a "shape theorem" for the red zone: lZ(t)/t 
converges with probability 1 to a deterministic shape. This proves the existence 
of a (maybe non isotropic) propagation velocity of the rumor or the combustion 
front. In this context C\ and C2 are respectively uniform lower and upper 
bounds of this possibly non isotropic front propagation velocity. 

v) It is believed that in the general case Dr > this propagation velocity does 
not depend on Db (see [S] - note 38). 

In this paper we give an upper bound on the propagation velocity, i.e., a 
ballistic upper bound on TZ(t) of the kind (|1.2|) with explicit dependence of 
C2 on the density p and no dependence on Db- This bound will be, in all 
dimensions, of order max(p, y/p). We argue that for d > 2 or p > 1 this bound 
can be optimal (in p), while for d = 1 and p < 1, we give in the simplest case 
Db = a better bound of order p. In addition we prove that our upper bound in 
max(p, y/p) holds for a larger class of models. We prove it, on the one hand, for 
those models in which red particles perform independent random walks while 
blue particles follow any kind of independent bistochastic process (see below). 
On the other hand, we give an analogous upper bound for models in which the 
rumor diffuses through a "contact process" inside an interacting particle system 
with exclusion and possible attraction (simple exclusion, Kawasaki dynamics) 
when a low density limit allows for a Quasi Random Walk approximation as 
introduced in [7j. 

1.2 One upper bound for many models 

We now define the first class of models we will work with. Like previously we 
start with a density p > of particles putting independently in each site z G Z d 
a Poissonian number of particles with mean p. We then put labels 1, 2, 3, ... 
on particles, we call Zi the position of the particle i and for alii > we will call 
Xi(t) € Z d and Yi(t) £ {R; B} its position and its color at time t. With each i 
we associate Zf and Zf two continuous time Markov processes on Z d in such 
a way that: 

• all these processes start from and are independent between them; 

• Zf is a simple random walk process with diffusion constant or jump rate 1; 

• Zf is a bistochastic process, i.e., satisfies 

Vz e Z d ,Vt > 0, p ( z o + Zf(t) =z) = l (1.3) 

This includes simple random walks with constant or time dependent jumps 
rates, long range random walks with drift, various deterministic pro- 
cesses, . . . 

• The Zf 's (like the Zf's) have the same law. 
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At time t = we choose one particle io with some probabilistic or determin- 
istic rule, we change the origin to put it where io stands, we give the red color 
to the particles in the new origin and the blue color to the other particles so 
that, for all i, 

Xi(0) = Zl -z lQ (1.4) 
Yi(0) = R ifX i (0) = (1.5) 
Yi(Q) = B if X t (0)^0 (1.6) 

Then each particle i follows the moves of Zf while Yi = B, turns red when 
it meets a red particle and then follows the moves of Zf-. More formally, with 
for all i, 



' if Xi(0) = 

inf{i > : Yi(t-) = B,3j ? i,^-(t_) = R,X t (t) = Xj(t)} 
if Xi(0) ^0 

with the usual convention inf = +oo, we have 



(1.7) 



XiW ~ \ X i (0)+ZB( Ti )+ZX(t-T i )tft>T i ^ 

We will call process of type RB any process that can be built in this way. 
The Kesten and Sidoravicius reaction-diffusion model is a process of type RB 

when Dr = 1. We will call it KS process. The general case Dr > can be 
mapped on the KS process by a simple time rescaling. 

Setting, like previously, for all f > 0, 

K(t) :={zeZ d : 3i > 1,3s e [0; t], {X tl Y,)(s) = {z,R)} (1.10) 
we will prove 

Theorem 1 There is a positive constant Sd that depends only on d and such 
that, for any RB process and for all t > 

p( 3 , eR(() \ B („,^£^L))<_ {1 . u) 

As a consequence, using the Borel-Cantelli lemma we get: 

Corollary 1.2.1 There is a positive constant 5d that depends only on d and 
such that for any RB process, with probability 1 

1l(t)ciB 0, Kt ' WH ' (1.12) 



5 d 

will hold for all t larger than some finite random time T 
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We will give an analogous result for a second class of models. In dimension 
d = 2 we consider a low density lattice gas, with density p, that evolves with 
exclusion and attraction inside a large finite box A(p) with periodic boundary 
conditions and according to the following Kawasaki dynamics at inverse tem- 
perature (3 > 0. With 

N:=p\A(p)\ (1.13) 

where |A(p)| denotes the volume of A(p), we will write f)i(t) S A(p) for the 
position at time t of the particle i in {1; . . . ; N} and rjt € {0; l} A ( p ) for the 
configuration of the occupied sites in A(p), in such a way that, for all t > 0, 

E Vt(z) = N (1.14) 

zSA(p) 

The energy of a configuration r\ € {0; l} A ( p ) is 

H(ri)--= J2 ~ Ur l(xMy) (1-15) 

{a;:iy}GA(p) 
|x-y|=l 

where | • | stands now for the Euclidean norm and — U < is the binding energy. 
With each particle we associate a Poissonian clock of intensity 1 . At each time 
t when a particle's clock rings we extract with uniform probability a nearest 
neighbor site of the particle, say i. If this site is occupied by another particle 
then i does not move. If not, we consider the configuration rf obtained by 
moving i to the vacant site, then with probability 

p = e -mw)-H(v)] + (L16) 

i moves to the vacant site and, with probability 1 — p, i remains where it was 
at time t_. Observe that the case U = corresponds to the simple exclusion 
process. 

In addition we choose at time t — some particle i$ according to some 
probabilistic or deterministic rule and give to i®, as well as to the particles that 
share with iq the same cluster at time t = 0, the red color, while all the other 
particles receive the blue color. Like previously a red particle will definitively 
remain red and a blue particle turns red as soon as it shares some cluster with 
some red particle. We call RBK process this dynamics and, for all t > 0, the 
red zone lZ(t) is defined like above. 

To control the propagation of the red particles in the regime p —> we will 
use the low density to reduce the problem to simple random walks estimates. 
This is more challenging when p and f3 go jointly to and +oo: in this case we 
have not only a low density regime but also a strong interaction regime. We 
will then deal with this more challenging regime only, setting p — e~ A/3 for A a 
positive parameter and sending (3 to infinity. We will write A^ for A(p) and we 
will choose |A^| = e ei3 for some real parameter > A. This regime was studied 
in [7] where a "Quasi Random Walk (QRW) property" was proved up to the 
first time of "anomalous concentration" T a ^\. For a a positive parameter that 
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can be chosen as close as as we want and A a slowly increasing and unbounded 
function such that 

A(/?)lnA(/3) =o(ln/3) (1.17) 

(for example A(/3) = ^/\n (3), T a \ is defined as the first time when there is a 
square box A C with volume less than eP^~ a l^ that contains more than 
A/4 particles. We will recall and use this QRW property to prove 

Theorem 2 For the RBK process, for all 5 > and all C > 0, uniformly in 
the starting configuration, and uniformly in T = T((3) < e c @ 

P (3z e 1Z(T) \ B (0, e 5l3 ^/pT) and T a x > T) < p~ 3 e 50 exp {-e~ 5f3 pt} + SES 

(1-18) 

where SES stands for "super exponentially small", i.e., for a positive function 
f that does not depend on T and the starting configuration and such that 

lim -ln/(^) = -oo (1.19) 

As a straightforward consequence: 

Corollary 1.2.2 For the RBK process, for all 5 > and all C2 > C\ > A, 

uniformly in the starting configuration, and uniformly in T — T((3) such that 
T > e Clfj and T < e C2f3 , 

P (3z € 1Z(T) \ B (0, e Sl3 y/p~T) and T QjA > T) < SES (1.20) 

Of course these results would be of no use if we were not able to have 
some control on T a \. But in [7] we discussed the fact starting from a "good 
configuration" T a ^\ is "very long". For example we proved that in the case 
A > 2U , starting from the canonical Gibbs measure associated with H, for all 
C > 0, 

P (T Q , A < e c/3 ) = SES (1.21) 

As a consequence of these results we will prove a long range decorrelation 
of dynamical events in this low density regime. Given A^ 1 ) and A^ 2 ) two square 
boxes contained in A^ we will denote by d(A^\ A^ 2 ') their Euclidean distance 
and by (J r f 1 '^)t>o an d {Ft^)t>o the filtrations generated by (j) tA -r a A | A (i> ) t > and 
(VtAT a Ja(2))*>o- With these notations: 

Theorem 3 For the Kawasaki dynamics, for all 5 > and all C > 0, uniformly 
in the starting configuration, uniformly in T — T(/3) < e c ^ , uniformly in AS 1 ' 
and A^ 2 ) such that 

rf(AW,A^) > e Sf> max(v / T, y/pT) (1.22) 
and uniformly in (A^\A^) G x 



P 



f A d) n a< 2 )) - p p (a< 2 )) I < SES (1.23) 

In the study of the low temperature metastable Kawasaki dynamics (the 
case U < A < 2U, see PQ) we will need such a long range decorrelation property 
(see [7])- This constituted the initial motivation of this paper. 
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1.3 How good are our bounds? 

In this paper we will not give any lower bound on the propagation velocity. But 
we give here some heuristic that indicates that max(p, ^fp) should be the right 
order of the velocity propagation in different situations. 

Consider for now the KS process in dimension d = 2 with p < 1 and in 
the special case Db = Dr = 1. lZ(t) should then look like a kind of ball 
that contains all the red particles and very few blue particles. In addition 
Db = Dr implies that, except for the color propagation, the particle system 
starts and remains at equilibrium. Let us call n(t) the number of red particles 
at time t. Since only the particles at the border of TZ(t) should contribute to 
the propagation of the rumor and since a particle typically waits for a time 1/p 
before meeting another particle, we should have 

dn ~ est \fnpdt (1-24) 

where 'est ' stands for a positive constant the value of which can change from 
line to line. As a consequence 

V^~cstpt (1.25) 

If r(t) stand for the radius of the smallest Euclidean ball that contains lZ(t) we 
should have 

n~cstr 2 p (1.26) 



so that 



r ~ est — ~ est ^/pt (1.27) 
\[f> 



If p > 1 we will typically have p particles per site and (|1.24ll turns into 



dn ~ est p. / — pdt (1-28) 
V P 

so that 

r ~ est ~ est pt (1.29) 

If d > 3 or Dr ^ Db we do not have such kind of heuristic. In the former 
case indeed TZ(i) should be a more complex fractal object, in the latter case the 
system does not stay at equilibrium. However Theorem Q] says that an upper 
bound of order max(/o, ^fp) holds independently of Db and the dimension. 

For d=l, Dr = Db and p < 1 the previous heuristic has to be modified. In 
this case the typical inter-particle distance is 1/p and a particle typically waits 
for a time 1/p 2 before meeting another particle. Then (|1.24jl and (|1.26ll turn 
into 

dn ~ est p 2 eft (1.30) 
n ~ est rp (1-31) 
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We will prove an upper bound of order p for the simplest case of the KS process, 
that is Db — 0, also known as frog model. 

Proposition 1.3.1 For the KS process in dimension 1 and with Db = 0, there 
is a positive constant 8 such that for all t > 



As previously we then get with the Borel-Cantelli lemma: 

Corollary 1.3.2 For the KS process in dimension 1 and with Db = there is 
a positive constant 5 such that, with probability 1 



will hold for all t larger than some finite random time T . 

We will give in section 2] some indications on how one can extend the simple 
proof of Proposition ll.3.ll to the general case of the KS processes. This is rather 
technical and we will not go beyond these indications. 

1.4 Notation and outline of the paper 

We will write est for a finite and positive constant that depends only on the 
dimension d and the value of which can change from line to line. Given d > 1 
we will write | • | for the d-dimensional Euclidean norm. Given a Markov process 
X and x in its state space, we will write P x for the law of the process that starts 
from x. 

In section[2]we prove simple random walk and large deviations estimates and 
we recall some definitions and properties regarding the QRW approximation for 
the Kawasaki dynamics. In section [3] we prove Theorem [T] for the frog model as 
well as Proposition 1 1 . 3 . Il In section 0] we prove Theorem [1] in the general case 
as well as Theorems [2] and [3] 




(1.33) 




(1.34) 
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2 Preliminaries 

2.1 Random walk and large deviation estimates 

Lemma 2.1.1 Let N and N' be two independent Poisson variables and 7 > 1 
such that E[N'] > jE[N]. Then 

i) P(N>jE[N]) <exp{-£[AT](7ln7- (7- 1))} (2.1) 
ii) P (n < ffl) < cxp {-E[N] ((l - i) - ^) } (2.2) 

/TV TV' \ 

" z) p {EW}- 1 W r u - 2exp{ ~ E[Nmnt ~ {t ~ m (2 ' 3) 

7 - 1 

wi/i i := — sll;7[ 

In 7 

Proof: We just use the Chebyshev exponential inequality. With A = E[N] we 
have, for any t > 0, 

P(N > 7A) < e^ x E[e tN ] = exp{-A(t 7 - (e 4 - 1))} (2.4) 

Optimizing in t we find (|2.ip with t = In 7. Similarly, for any t > 0, 

P(iV < A/7) < e txh E[e~ tN ] = exp{-A((l - e"*) - t/7)} (2-5) 

Optimizing in t we find (|2.2p with t = In 7. Finally we have, for any t > 0, 

^ > T^j) < P (N > + P (V < 1 W ) (2.6) 

By (j2~Tj) and (pE?)) this gives, if < > 1 and t < 7, 

TV AT' 
>7" 



E[N] ~ ' E[N'[ 
< exp{-A(tlni- (t- 1))} 

+ exp{-A7((l-^)+^ln^n (2-7) 



The two terms of this sum are equal when 



7-1 

t=-\ (2.8) 

ln7 

The concavity of the logarithm ensures 

1 < -In- =1117 < 7- 1 (2.9) 

7 7 

so that 1 < t < 7 and this gives (|2.3[) . □ 
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Lemma 2.1.2 Let C, be a d- dimensional continuous time simple random walk 
with rate jump 1. For all t > and z (zlf 1 



• */ \z\ < t then 



if \z\ > t then 



Pom 



. est 



cst|z| 2 



P (((t) = z) < est exp{-cst|2;|} 



(2.10) 



(2.11) 



Remark: Since we just need an upper bound on these probabilities we do not 
need the usual condition \z\ = o(t 2 / 3 ) of the local central limit theorem. How- 
ever, working with continuous time random walks, we have to treat separately 
the case \z\ > t. 

Proof of the lemma: We will prove slightly different but equivalent estimates: 
(pHI))) when \z\ < 2t (|2"TTj) when \z\ > 2t. 

For the case \z\ > It we apply the previous lemma. If £ reaches z in time 
t then the number of its clock rings up to time t is larger than or equal to \z\. 
Since this number has a Poissonian distribution of mean i, this occurs, by (|2.1[) . 
with a probability smaller than 



exp < —t 



N ln N 
t t 



t 



i 



< exp 



est | z I 



t 



-est \z\ 



(2.12) 



(for the last inequality we used that \z\jt was bounded away from 1.) 

For the case \z\ < 2t we first observe that, working with a continuous time 
process with independent coordinates, it is enough to prove the result for d = 
1. Then we prove the estimate for £ the discrete time version of such a one 
dimensional process. Without loss of generality we can assume that z G Z is 
non negative. If z < n/2, then, by the Stirling formula, 



Pi 



b (C(n) 



< 



< 



< 



est 



2" 



est 



' n—z ' 
v 2 / 



est 



exp {— nl(z/n)} 



1 - 



(2.13) 

(2.14) 
(2.15) 



with 



T, N 1 + * 1 /■- N 1 

:=^^ln(l + x) + - 



•ln(l-ar), xe[-l;l] 



It is immediate to check that 



J(0) = J'(0) = 
V ! c€]-l;l[,/ // (x) = I ^ P r>l 



(2.16) 



(2.17) 
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As a consequence, for all x £ [— 1; 1], 



m > y 

and this gives, for z < n/2, 

Po(((n) = z) < ^lexp- 
This is easily extended to the case z > n/2, i.e., z/n > 1/2: 
Pa K( n ) — zj < cs t exp {— nl(z/n)} 



z 

2n 



(2.18) 

(2.19) 

(2.20) 
(2.21) 

(2.22) 

(2.23) 

(2.24) 



where JV is a Poissonian variable of mean t. By (|2.ip . (|2.2p applied with a large 
enough 7 we can find two positive constants c±, C2 with 4ci < C2 such that 



< est exp <^ -n- 8/(1/2) — - 
' 2n z 



< cst A / — exp<j - — 



< 



est 



6XP ^2^ 



Finally we use the previous lemma to prove (|2.10[) . We have 
Po(({n)=z) < E 



est 




— = exp < 

_Vn I 


. 27V }_ 



est 



Pq (C( n ) = z) < —jr exp \ -ci— \ + exp {-2c 2 t} 
est 



< —j= ( exp \ -c\ — \ + exp {-C2O 



and we get (|2.10p using z < 2t, i.e., At > z 2 /t. 



(2.25) 
(2.26) 
□ 



2.2 Quasi Random Walks 

With the notation we introduced in section 11.21 for the Kawasaki dynamics and 
given an arbitrarily small parameter a > as well as an unbounded slowly 
increasing function A satisfying (| 1 . 1 7|) . we recall in this section a few definitions 
and results from [7|. 

Definition 2.2.1 A process Z = [Z\\ . . . ; Z^) on Ap is called a random walk 
with pauses (RWP) associated with the stopping times 

= a i;0 = T i>0 < <j iA < n,i < < n t2 < ■ ■ ■ iE {!,..., N} (2.27) 
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if for any i in {1; . . . ; N}, Z$ is constant on all time intervals [<7i,fc, Tj ; fc], fc > 0, 
and i/ i/ie process Z = [Z\ , . . . , Zjv) obtained from Z; by cutting off these pauses 
intervals, i.e., with 

Zi(s) :=Ziis+ T i,k - <n,k , s > (2.28) 

V k<Ms) J 

where 

Ms) := inf J j > : s + ^ r a - a,, fc < o^- I (2.29) 

[ k<j J 

is an independent random walk process in law. 
Now with 

T a := e ( A - a ^ (2.30) 
Quasi Random Walk processes arc defined as follows. 

Definition 2.2.2 We say that a process £ = (£i,...,£jv) 071 is a Quasi 
Random Walk process with parameter a > up to a stopping time T, written 
QRW(a,T), if there exists a coupling between!; and a RWP process Z associated 
with stopping times 

= <7i,0 = Ufi < <7 i; l < Tis < C7 i;2 < Ti.2 < . . . i e {1; . . . ; N} (2.31) 

such that £(0) = Z(0), for any i in {1, ... , TV} & and evolves jointly — Zi 
is constant) outside the pause intervals [a^fc, r^fc], k >0, and for any to > the 
following events occur with probability 1 — SES uniformly in i and to: 

Fi(t ) := {»{fc>0: <T iik e [t Q A T,(t + T a ) AT}} < l(p)} (2.32) 

Gi(to) := {vfc>0,Vt>to,o-i, fc e[toAT,(to+T a )AT] 

ICC* A T »,fc A T )«(* A A T )l < (2-33) 

/or some f} i— » Z(/3) t/iai satisfies 

lim i lnZ(/3) = (2.34) 

In words, the fact that for each i the events Fj(fo) and Gi(t ) occur for all 
to > means, on the one hand, that in each time interval before time T and 
of length 1/ p almost, there are few pauses for the associated RWP Zi (a non 
exponentially large number) and, on the other hand, that £j stays close to Zi in 
the sense that during each of these few pause intervals the distance between the 
two processes cannot increase of more than the same non exponentially large 
quantity /. 
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Proposition 2.2.3 For any unbounded and slowly increasing function A that 
satisfies \1. 17\ ) and any positive a < A, rj is a QRW(a,T a \) process. 

We refer to [7J for the proof. As a consequence of this QRW property we 
have for all <5 > 0, uniformly in the initial configuration and uniformly in T = 

T((3)e [2,T Q 2 ] 

P {l~a,\ > T,3t (z [0, T],3i G {1, . . . , N}, \fji(t) — ?7i(0)| > e sp T) < SES (2.35) 

In [7J we also introduced at any time to > a partition of {1; ... ; N} in clouds 
of potentially interacting particles on time scale T a : we associate with each 
particle i a ball centered at its position at time to with radius 

r:=e^y/T^ (2.36) 

we call Bq their union 

B Q :=U t B(r} t (t Q ),r) (2.37) 

and we say that two particles are in the same cloud if there are, at time to, in 
the same connected component of Bo- It is easy to check that if to < T a \ then 
no cloud contains more than A particles. And, as a consequence of (|2.35[) . with 
probability 1 — SES interactions between particles during the time interval 
[to, (to + T a ) A T a ,\[ will only take place inside the different clouds (and not 
between particles of different clouds.) 

3 The frog model 

3.1 Proof of Theorem [j] for the KS process with D B = 

There is a natural notion of generation in the model. We say that the first 
particle in the origin is of first generation and that a particle that turns red when 
it encounters a particle of fcth generation is of (k + l)th generation. (If a blue 
particle moves on a site with more than one red particles then its generation 
number is determined by the lowest generation number of the red particles.) 
Now, to drive the red color outside a ball an Euclidean ball -6(0, r) in time t, 
the first particle initially in z\ = has to activate at some time ti and in some 
site Z2 a second generation particle, this particle has to activate at some time 
t\ + ti and in some site 2:3 a third generation particle,. . .and, for some n, an 
nth generation will have to reach some site outside B(0, r) at some time 

ti H h t n < t. Taking into account the fact that more than one blue particle 

can stand in a site reached by a red particle ad using Lemma 12.1.21 we get, for 
all r and t, 

P(3zeK(t),\z\>r)<Q(r,t) (3.1) 
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with 

QM) -EE/ < e n«-'£* 

„>1 Zi,...,z„ + i ■ / *l + -+tn<t j 2 ,...,j n >0fe=2 J 
21=0 

z„+igB(0,r) 

n 5 Vcste-^l^-^lJ dtfe (3.2) 

where here like in the sequel we did not write, to alleviate the notation, that 
the integral is restricted to positive variables only. 

Permuting the last sum with the product, making a spherical change of 
variable and using the triangular inequality we get 

Q{r,t)<Y, L+---+r n >rP n ~ i n<li{r k ,t k )Vq2{r k )r d k - l dr k dt k (3.3) 
n > 1 Jt 1 +-+t n <t k=1 

with 

. . est cst r l . 

q\(r k ,t k ) := -^e ** (3.4) 

q 2 (r k ,t k ) = q 2 (r k ) := c S te~ cstr « (3.5) 

Grouping together the different terms according to the respective values of q\ 
and q 2 and using 

' n 1 



< 2™ (3.6) 



we get 

Q(r,t) 

< 



\ L+R 2 >rJ2J2 (f) { L+-+r j >R 1 P i Y[<li( r k,tk)ri l dr k dt k \ 

PJT 1 +T 2 <t n > lj=0 \J/\Jti+-+t j <T 1 k=1 j 

( / 1+ ...+r„_ 3 > fe P n ~ j II 92(r k )ri- 1 dr k dt k ) dRidRadTidT 2 (3.7) 

\J t 1 +-+t n - j <T 3 k=1 J 

<- L+R 2 >rQi{Ri,T 1 )Q 2 (R2,T 2 )dR 1 dR 2 dT 1 dT 2 (3.8) 

P J T,+T 2 <t 



7i+T 2 <t 

with for j = 1, 2 



« n 

Qj{Rj,Tj) := r 1+ ... + r n >R j ^p) n llQj(rk,t k )r d k - 1 dr k dt k (3.9) 

n>l^*H l-*n<Tj fc=1 

For any i?, T > we will estimate separately <2i(.ffi,Ti) and Q 2 {R 2l T 2 ). 
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We have 

Q 1 (R,T)<J2(^P) n L... + r n > R f[e- cst ^ fj^Y''^ ( 3.io) 

„>! J ti+-+t„<T fe=1 Wife/ VEfe 

Making a change of variable = cstr^/^. and observing that, by the 
Cauchy-Schwartz inequality, 

E k VtkV^k> R ( J2 k x k > est R 2 /T 
we get, with T the Euler function, 



Ql (R,T) < EMri"L.. +I „>c«VTll^^^( 3 ' 12 ) 

J t-t-l kt„<T I — i Xi. 

_ 

r(d/2) 



n>l ti+-+t»<T fe=i 



Since the volume of the n-dimensional simplex of side-length T is T n jn\ and 
the sum of independent variables with a T distribution follows a T law, 

Qr(R,T) < e-^l--^ (3.14) 

4^ n! Jm>cMtB?/T ring) 



n>1 n - Jx>cstR 2 /T 

nd 

—<—F | TV' < 

n! 



(3.15) 

n>l 

< e cstpT P(7V' < est TV) (3.16) 



where AT and N' are independent Poissonian variables of mean est R 2 /T respec- 
tively. Now, for any large enough 7, if R > jyfpT, then by (|2.3|) 

W*1 s ^K^^'i™) (3 - 17) 

so that, for any large enough 7, 

Qi(fl,T) < e cst " T exp{-cstVpi?l [7V?Ti+oo[ (i?)} (3.20) 
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Turning to Q2(R,T) we have 



Q2(R,T) < £(cstp) n L + ... +rn > R l[e- cstr *r d k - 1 dr k dt k (3.21) 

n>l Jt 1 +-+t n <T k=1 

„ n 

^Mp)" L+...+x„>c B tflII e ~ 1B * a! fc~ lda: * dtfc ( 3 ' 22 ) 
n>l J h+-+t n <T k=1 



< 



< 



y^ETj e~V^-^- (3.23) 
^ »l J x >cstR r(nd) 



n>l 



< e cstpT P(N' < est iV) (3.24) 

where iV and iV' are independent Poissonian variables of mean est i? respectively. 
Then, for any large enough 7, if R > jpT, we get by (12. 3[) 



*M S e «- p (_«. (3,5, 

so that, for any large enough 7, 

Q 2 (i?,T) < e cstpT exp{-csti?l [7pT:+oo[ (i?)} (3.27) 
Turning back to Q(r,t), we get, for any large enough 7, 
Q(r,t) 

< I / cstp(Ti+T 2 ) 
P J Ti_+T 2 <t 

exp{-cst (y/p~RilL h ^p Tu+oc[ (Ri) + i? 2 l [7pT2j+oo[ (i? 2 ))} 

dR 1 dR 2 dT 1 dT 2 (3.28) 

<• i- /" est pt 

P J Tx+T 2 <t 

exp{-cst (y/pR\$.ft l a u+00 [(yfpR-i) +i? 2 !L[ 7 pT 2 ,+oo[(-R2))} 

dR 1 dR 2 dT 1 dT 2 (3.29) 

Now if p < 1, then 

^2H[ 7 pT 2 ,+oo[(i?2) > V / P-R2l[ 7 pT 2 ,+oo[(\/P^2) (3.30) 

and if p > 1 , then 

(y/pRl) > flii [7pTl , +oo[ (Ei) (3.31) 

As a consequence, with 

:= max(/9, -Vp) and = — R4, i = 1, 2 (3.32) 
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we have 

Q(r,t) 



< 1 / est pt 
"'" J T t +T 2 <t 

exp{-cst (X 1 1L hpTu+oo[ (X 1 ) + X 2 H hpT2i+oo[ (X 2 ))} 

dX 1 dX 2 dT 1 dT 2 (3.33) 

est /yt p 

< £ / „-cst(Xi+J>f2-7P(Ti+T2)) 

- ~T~ X 1 +X 2 >pr/p e 

' J Ti+T 2 <t 

dX 1 dX 2 dT 1 dT 2 (3.34) 



If r > 2jpt, i.e. 



£ > (3.35) 
2p 



then 



QM) < — — L+x^/fa-^-^dXidXidTidTi (3.36) 

" •/ Ti+T 2 <t 
„cst pt 

< cst ^ ^ e -c St f, (337) 

p p 

< est— ^e" cst ^ (3.38) 

P 



est 
7 

and, with a large enough 7, we get 

est 



Q(r,t) < r^ e — t P t (340) 



P 4 



□ 



3.2 Proof of Proposition 11.3.11 

In the previous proof we could have use, instead of the estimates from Lemma 
12. 1.21 on Po(((t) = z)dt, an estimate on 

dP (t,(C) < t) = P (t z (C) e [t,t + dt}) (3.41) 

with 

7>(C) :=inf{i>0: 0) = z} (3.42) 

While in dimension d > 2 the two quantities are quite close, in dimension d — 1 
they are substantially different. In addition, using r z (C) in dimension 1 allows 
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for a simpler proof of a stronger result when p is small enough. Indeed, for all 
r and t, 

P(3z€ll(t),\z\ >r) 

Y[dP a (r rk (C)<t k ) (3.43) 

n^i nt...r»dr ~ lH ^*"<* fc=l 

^ E' n_1 E E fl>wo<<) (3-44) 

n>l i?>r ri + ...r„=fl 

Then, by the reflexion principle and Lemma 12.1.21 

est x x , p n B!> 



P(3ze1l(t),\z\>r) < ^J2J2^-( e ~ CStR2/Tve ~ CStR )^ 

P R>rn>l ' 

-J2e pR (e- cstR2/T ye- catR ) (3.46) 



est 

< 



R>r 

Now if r > jpt for some large enough 7 we get, for p small enough, 

est 

P(3zeK(t),\z\>r) < J2e- cstpR (3.47) 

P R>r 



CS ^ -est pr 
P 2 ' 



(3.48) 



< ^e-^" 2 ' (3.49) 
P 

This proves Proposition 1 1 . 3 ,T1 for small p's. When p is bounded away from 0, 
Proposition II. 3. ll is just a consequence of Theorem [1] for the frog model. 

4 RB and RBK processes 
4.1 Proof of theorem [I] 

We can proceed like in the case of the frog model except for the fact that a 
particle does not anymore turns red at the same point where it started. We 
have then to sum on the possible starting points. With the notation 

s fe =*! + •• ■ + t k - 1 , k>2 (4.1) 

and for any i > 1 we have 

P(3z€ll(t),\z\>r) (4.2) 

<E E / E f[e- p QhP(4 + Z*( Sk )=z k ) 

z 2 ,...,z n 



n>l zi,...,z»+i "«i z i z i k=2 

2:1=0 

,+ i£B(0,r) 



Zl ~° 32,...,i»>0 



/ CQ"t~ est — z h,\^ \ 

II ^7J e 5 Vcste" 08 * 1 ^ 1 -^ 1 dt k (4.3) 



fc=i \*& 
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Now permuting the last sum with the product and using (|1.3[) we get 

P (3z G n(t),\z\ > r) < Q(r,t) (4.4) 

with Q(r,t) defined in (|3.2[) and estimated in the previous section. □ 
Remark: Unfortunately the proof of Proposition 11.3.11 cannot be extended so 
simply to the general case, even if we restrict ourselves to KS processes. To do 
so we would have to link the differential 

dP (t Zr (Cr) <t)=P (t Zr {(r) e [t,t + dt\) (4.5) 

with the sum 

P (o^-b) MO? - Cr) e [t, t + dt], Ch(«) = zr) (4.6) 

Zfl>0 

with £r and Cb independent continuous time random walks with jump rates 
Dr = 1 and > 0. In the case Db = 1 this can be done using the indepen- 
dence between Cb — Cr an d Cb + Cr- I n the case -Db ^ 1 we can only use an 
"asymptotic independence" between Cb — Cr an d Cb +DbCr- I n both cases this 
is a quite technical task: we will not go in this paper beyond the result for the 
frog model. 



4.2 Proof of Theorem [2] 

We can adapt the proof for the frog model using the QRW property and the 
last observations of section [ 



P (3z g K(T),\z\ >i?,T>T Q , A ) 

\ \IT /T a ~\ 

n=l z„ i1 fA,/ tl+ " ,+t »i i fcl 



< 



zi,...,z„ + ieA fj 

21=0 



3;2 




In this formula the first sum is limited to \\IT /T a ~\ since T is at most exponential 
in [3 and in each interval of length T a , with probability 1 — SES, interactions 
are limited to clouds that contains A particles at most and particles are coupled 
with random walks with I pauses at most. The factor I 2 is due to the fact that, 
with probability 1 — SES, in each pause interval the distance between a particle 
and its associated random walk with pauses increases of I at most, one factor A 
is due to the fact that A red particles at most can leave a given cluster before 
T a \ and the last factor A 2 is due to the fact that at each time t < T a \ a given 
particle can turn red other particles inside a radius A at most. 

Then we can repeat the calculation of section [3~T1 with two main differences. 
On the one hand we do not have anymore the factor p™ -1 in our sum, on the 
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other hand this sum is limited to \XIT/T a ~\ . Instead of (|2.3p we use then 
repeatedly. For example defining Qi and Qi in an analogous way and observing 
that for any 6 > 0, A and I are smaller e s @ for j3 large enough, we have now 

Qi(R,T) < V ^ — r-tP (N' < e 5/3 P T) + SES (4.8) 
^ — ' nl 

n=l 

with N' a Poisson variable of mean est R 2 /T. For any <5 X > 5/2, if R > e Sl0 \fpT 
the last probability can be estimated from above by 

P (N' < e sp pT) < exp {-est pt} + SES (4.9) 

while the last sum can be estimated from above by 

r-^p^i (e sp T) 

V ^ — < exp{e^T}P (N < est e S0 pT) + SES (4.10) 

n=l 



with N a Poisson variable of mean e ' pT, so that, by (12. 2p . 

V 1 < exp{e 2 ^pT} + S£S (4.11) 



n! 

n— 1 



Putting everything together we get, for any R, T, 

Qi(ii,T) < eMe 2 ^pT}e X p{-cst^pR1L [e!l ,^ +oo[ (R)} + SES (4.12) 

We can estimate Q2 in the same way and the rest of the calculation goes like in 
section 13.11 □ 



4.3 Proof of Theorem U 

Given A (1) and A (2) with 

d(A«,AW) > e 5/3 max(Vr,Vpr) (4.13) 
we define a new coloring process. With 

5:=A«UA( 2 ) (4.14) 

and 

W:=i^zeAp: inf |z - b\ > e - 5l3/2 d(A^\ A< 2 >)J (4.15) 

we say that all the particles that start from B are black, all the particles that 
start from W are white and all the particles that start from (B U W) c do not 
have any color at time t — Then, for t > 0, black particles keep their black 
color, white particles keep their white color, non-colored particles that enter 
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B turn black, non-colored particles that enter W turn white, and non-colored 
particles that share some cluster with a colored particle turn black or white 
choosing randomly a colored particle inside the cluster and taking the same 
color. We can define a black zone and a white zone like we defined the red zone. 
As a consequence of corollary 11.2.21 with probability 1 — SES, the black and 
white zones will not intersect up to time T A Ta.,\ and we will never see black 
and white particles in a same cluster up to time T A T a ,\- 

Now we couple in the more natural way the previous process, with a process 
that starts from the same initial configuration, uses the same marks and clocks 
for the particles and evolves in the same way except for the fact that each particle 
in W or that enters in W disappears. For this process the restrictions of the 
dynamics to A^ 1 ) and A^ 2 ) are clearly independent and the previous observation 
shows that, with probability 1 — SES, these restrictions for the two processes 
coincide up to time T A T a ,\- This proves the theorem. □ 
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